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ABSTRACT 

The generator QEDPS, a parton shower model in QED, developed by the present 
authors has been improved. By a careful study on the shower algorithm it was found 
that some finite contributions have been missed in the original version. They are 
small, but cannot be neglected when the contribution of the soft photons is domi- 
nant. A method is presented to correct these finite pieces and then a new generator, 
improved QEDPS, is proposed, which is able to make more precise prediction for the 
processes observed in e + e~ annihilation as far as the initial radiation is concerned. 
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Section 1. Introduction 

In detailed studies of radiative corrections to e + e~ processes one needs Monte 
Carlo generators. Several authors have proposed such generators so far JT]] — 0. We 
have already published QEDPS for radiative corrections in the e + e~ annihilation |2| 
and also in the Bhabha scattering ||. These generators are a QED version of the 
parton shower which has been originally developed for the study of QCDJ3J]. The 
algorithm is based on the Q 2 -evolution according to the renormalization group equa- 
tions, which enable us to sum up all the collinear singularities. 

The current version of QEDPS is written in the leading-logarithm(LL) approx- 
imation. Recent experiments at e + e~ colliders, however, can provide very precise 
data and this in turn demands that the theoretical uncertainty in a generator should 
be less than 1% level. In order to respond to this requirement we have been study- 
ing the complete next-to-leading-logarithm of the QEDPS like QCD||. During this 
study we found that some finite contributions had been neglected in the original 
version. This takes place in the algorithm we adopted. These contributions are 
small, but actually cannot be ignored for the observables for which soft photons are 
important. 

We explain how these contributions arises and how to do with them in the next 
section. Their origin is directly connected to the kinematics used in the generation 
algorithm. First we take the single cascade scheme which means that the electron 
with the momentum parallel to the axial gauge vector does not make showers. In 
section 3 the same problem will be discussed for the double cascade scheme. In this 
case any electron is able to make showers. Then one possible way to correct them 
is proposed and some results by the improved generator will be presented in section 
4. Final section is devoted to summary and discussion. 

Section 2. Q 2 -evolution 

The algorithm of the parton showers is based on the Altarelli-Parisi(AP) equation 
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which governs the Q 2 -evolution of the structure function of quarks and gluons in 
QCD. In QED the strong coupling a s is replaced by that of QED a and the color fac- 
tors are removed from the AP equations. The equation for non-singlet (^-evolution 
is given by 

dD(x,Q 2 ) a [ l dy 2 

Here D(x,Q 2 ) is the structure function of the electrons with x being the electron 
momentum fraction and P+(x) = ((1 + x 2 )/(l — x))+ is the split function. In the 
following we shall neglect the running effect of a only for the sake of simplicity. By 
taking the moments of this equation, one can easily obtain the solution 

D(n, Q 2 ) = jf 1 dxx n - 1 D(x, Q 2 ) = A n exp (~^7n log(Q 2 /Qo)) , (2) 

where A n is the integration constant and 

7 „./ o ^P +W .----_-2g J . (3) 

is the anomalous dimension. Thus the Q 2 -evolution of the structure function is 
determined by 

D(n, Ql)/D(n, Q 2 2 ) = exp (|- 7 „ log(g 2 /Q 2 )) . (4) 

Let us show the algorithm how to generate the parton showers. First we have 
to regard Eq.(fJ) as an equation which describes the stochastic process of emitting 
photons. For this we modify the split function as follows 0; 

P+(x) =6(l-e- x)P(x) - 5(1 - x) f 1 6 dyP(x), P(x) = i±£L (5) 

Jo 1 — X 

Here e is a small quantity and 9 is the step function. Then the evolution equation 
can be cast into the form 

D(x,Q 2 ) = U(Q 2 ,Q 2 s )D(x,Q 2 s ) 

+ f /T ^ n ^ 2 ' -ny)D{x/y, K% (6) 
27r Jq 2 s K J x y 
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where Q 2 is the minimum of Q 2 and U(Q 2 ,Q' 2 ) is the non-emission probability 
defined by 

W,rt = °*(-££%C «(•))• P) 

The evolution equation in the integral form Eq.(|D allows one to take it as that for 
stochastic process and the algorithm of the photon shower consists of the following 
steps. 

(1) Set Xb — 1, where Xb becomes the fraction of the light-cone momentum of the 
electron after the end of the shower or just before the annihilation. 

(2) If a given random number rj is smaller than H(Q 2 , Q 2 ), the evolution stops. If 
not, one finds the virtuality K 2 that satisfies rj = II (K 2 , Q 2 ) with which a branching 
is made. 

(3) Fix x according to the probability P(x) between and 1 — e. Then Xb is 
replaced by XbX. One should go to (2) by substituting K 2 into Q 2 and repeat until 
it stops. 

According to this algorithm the probability that no photons is emitted during 
the evolution from Q 2 to Q 2 is U(Q 2 ,Ql). The probability for the single photon 
emission is 

rQ 2 n 1 rl-e 

/ dK 2 U(K 2 ,Q 2 )^— / dxP(x)U(Q 2 ,K 2 ) 
JQl 2tc K jo 

rQ 2 a 1 r l ~ e 
= n(Q 2 ,Q*) j Q2 dK 2 -— j o dxP(x), (8) 

and for the double photons it is 

rQ 2 r K i nt 1 r l ~ e 

/ dK 2 2 d Klll{KlQ 2 s )^— 2 dx 1 P(x 1 )U(KlK 2 l ) 

JQs Jq 2 s 2ir Kf Jo 

x 7^172 / dx 2 P{x 2 )Tl{Q 2 ,Kl) 
Ztt K 2 Jo 

= U(Q 2 ,Q 2 ) / dK 2 — 2 / dx 2 P(x 2 ) / dK 2 ® — 2 / dx 1 P(x 1 ). 
JQl 27r A| Jo JQ 2 S Zt: Kf Jo 

(9) 

For the emission of any number of photons, we have similar expressions. 
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When one wants to see the distributions on x&, the (^-functions are inserted into 
the integrands, 

8(x b — x^.-Xn). (10) 
If one takes the moments, one gets a compact expression. 

D(n,Q 2 )=exp(-j Q2 — jf dxP{x){^ x - 1)J , (11) 
by noting that 

/ dx b x r ^ 1 5(x b - xix 2 ■ ■ -x N ) = (x 1 x 2 ■ ■ -xat)™" 1 , (12) 
Jo 

pQ 2 pK 2 pK 2 pK 2 1 N fQ 2 

L iKl L dK "--- L iKl L iKl= ^L dKl (13) 

If we assume a very small constant value eo for the parameter e, we reproduce 
Eqs.flU) because the corresponding structure function becomes 

D (n,Q 2 ) = e W (j-\og(Q 2 /Q 2 s ) J^ e ° dxPix^ 1 (14) 
= exp(|-log(Q 2 /Q^ + 0(6o)). (15) 

No problem arises, if eo is a constant. In the parton shower algorithm to generate 
events, however, we have _fT 2 -dependent e(see ref.0). In QEDPS it is given by 

e K = Ql/K 2 , (16) 

which comes from the fact that the kinematical boundary restricts the fraction x of 
the light-cone variable. Here Ql is a cutoff, a fictitious mass of the emitted photon 
introduced to regulate the infrared divergence. 

In order to see whether the Monte Carlo simulation with the above mentioned 
algorithm can reproduce the moment given by an analytic expression Eq.(^) we 
make a numerical comparison. We take Ql = Q 2 = 10 _6 GeV 2 ,Q 2 = 10 2 GeV 2 
and a = 1/20. The last value is assumed only to magnify the finite contributions. 
Figure 1 shows the both results from the simulation(points with diamond) together 
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with those from the analytic formula(solid curve). Apparently there is a sizeable 
discrepancy However this does not imply that the cutoff e# is wrong. 

This can be understood by looking at the ratio of the moments D(n, Q 2 )/D(n, Qf). 
This agrees completely with the analytic results. Hence one concludes that the struc- 
ture function with the cutoff = Q\jK 2 should give some finite contributions. As 
long as the Q 2 -evolution(the ratio at two different Q 2, s) is concerned, everything is 
correct. However if one is interested in the absolute value of the structure function, 
the algorithm is not sufficiently accurate and should be improved. The generated 
D(n,Q 2 ) must be able to reproduce the results of perturbative calculation which 
have been known in the literature as mentioned in Section 4. Thus one can see 
that in the actual applications the finite contribution does not cause any problem 
in QCD, but in QED more careful study is required. 

Let us now evaluate the finite contribution. First we define a structure function 
D s (n, Q 2 ) with the cutoff as 

/ n rQ 2 dK 2 r l ~^K \ 
D s (n,Q 2 ) ee exp ( — — jf dxP(x)(x n ^ - 1)J , (17) 

(Here we have attached a subscript s to the structure function to stress that it is 
defined in the single cascade scheme whose meaning will become clear in the next 
section.) When Ql <C Q 2 , we can modify 

D s (n,Q 2 ) = ex V J^j J\zP(1 - z)[(l - zT- 1 - 1)]^ 

= e W (^j^dz\og{^)P{l-z)[{l-z) n - l -l)^ (18) 

with 6q = Ql/Q 2 . From the last expression we find 

D s (n,Q 2 ) = D(n,Q 2 )D sf (n) + 0(e Q ), 

D sf (n) = exp(j^£dxP(x)\og(l-x)(x n - l -l)), (19) 

where D(n, Q 2 ) is given by Eq.(p|) with e = and D s f(n) is the finite contribution 
we have looked for. In x-space, the equation ([19]) becomes a convolution integral 
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with respect to x as 



D s (x,Q 2 ) = D{x,Q 2 )®D sf {x). 



(20) 



In figure 1 one can see that D s (n,Q 2 ) is completely reproduced by the Monte 
Carlo simulation when Q 2 S = Qq. If Ql <C Q 2 <C Q 2 , however, the finite contributions 
vanish. This implies that D(n,Q 2 ) should be coincides with the simulation(points 
with circle). This is also demonstrated in the figure. 

Section 3. The double cascade scheme 

The algorithm discussed in the previous section is for the case of the single cascade 
scheme in e + e~ annihilation. The characteristic of this scheme is that the vector 
to define the axial gauge is set parallel to the positron and as a result the shower 
develops only on the electron while the positron stays inactive. 

In the program for the actual generator, however, we adopt the double cascade 
scheme 0. This is another formulation which is more suitable than the single cascade 
in order to assure the symmetric treatment of the radiation from e + and e~ as both 
of these are able to develop showers. In this section we discuss this scheme in some 
details and look at what finite contribution comes out. 

First we give a review on the double cascade scheme, where photons are radiated 
parallel to its parent electron or positron but not anti-parallel to it. To make the 
argument clear, we take the simple case that only a single soft photon is emitted. 
The momentum of the photon is denoted as k and those of two (on-shell) electrons 
are Pi and P%, respectively. For simplicity we neglect the electron mass, that is 
P 2 = P| = 0. The soft photon contribution / is then given by 




(21) 



where (ab) = ao&o — a • b. We use the light-cone momentum defined as 




(22) 
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and choose a frame in which we have 

P lM =(P 1+ ,0_,O T ), P 2M = (0 + ,P 2 _,0 T ), 2(P l P 2 ) = Q 2 . (23) 

Now two new variables z and t, the fractions of the light-cone momentum of elec- 
trons, are introduced by 

_ _ (fcP 2 ) _ k + t= (fcPQ = (24) 



(PiP 2 ) Pi+' (PiP 2 ) P 2 _ 

In terms of these variables the soft photon contribution I is rewritten as 

J = ^/ 1 ^/ 1 i (25) 

47T Jo Z Jo t 

This integral is divergent so that a cutoff is needed. For this we impose the lower 
limit to the transverse momentum squared of the photon, Qq. 

k 2 T = z2(kP 1 ) = ztQ 2 > Q 2 Q , 

I{Ql) = ^f^f^-W). (26) 
An Jo z Jo t 

The events can be generated using these z, t in the full phase space. This generation 
method is the single cascade scheme. 

In the CM. frame we have P 1+ = P 2 _ = y / Q 2 /2 so that the z-component of the 
photon momentum is given by 

kz = K^ = ^L (z _ t) (27) 

This shows that the photon with z > t is emitted parallel to the electron, while 
that with t > z parallel to the positron. Consequently it would be understood that 
we can select the parent of a shower by looking at which of z or t is greater than 
the other. Thus a shower can develop from any electron independently and emitted 
photons are parallel to the electron or positron from which it has branched if one 
imposes the restriction z > t. This is the double cascade scheme. 

The parton showers are based on the renormalization group equations by which 
some physical quantity is calculated. It depends on the process considered and for 
the case of the deep inelastic scattering, it is the distribution on the Bjorken variable. 
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The parton showers must give the same distribution which can be calculated in the 
analytic approach. 

In the single cascade scheme the fraction of the light-cone variables is equal to 
the Bjorken variable. In the double cascade scheme this is not valid and we must 
find what combination of which variables generated in parton showers is equal to 
the Bjorken variable. 

Using the example of the single soft photon emission, we will look at the problem 
in details. Then the variable x B j can be introduced by the following way: 

Pi + q = P 2 + k, 2(P ig )+q 2 = 2(P 2 k), (28) 

J_ = 2(fld = 1+ Cft*). 
x Bj -q 2 -q 2 

Since -q 2 = 2(P 1 P 2 )(1 +t-z),we find 

x B j = l-z/(l + t) = l-z + 0(zt). (30) 

For small t the Bjorken variable is then determined by only z. This kinematics 
corresponds to the single cascade scheme, 

D(x Bj ) = 8(x Bj -l)(l-I(Q 2 )) 

a ri dz ri dt _ Q 2 /Q 2 )5{ _ 1 + z y 

47T JO Z JO t 

(31) 

In the double cascade scheme, on the other hand, the inequality between z and t 
directly connected to the choice that either of which electrons emits the photon. In 
other words, when z < t or (P 2 k) < (P±k), the role of z and t should be interchanged. 
Hence in our example of the single soft photon emission, the distribution over x B j 
becomes 

D(x Bj ) = S(x Bj -l)(l-I(Q 2 )) 

a ri dz, /-i dt, _ Ql/ Q 2 )e{zi _ tl)S{ -i + Zl ) 

47T JO Z\ JO t\ 

a /-i dz, ri dt2 0(z2t2 _ Ql/Q 2 )e{z2 _ h)5{x -l + t 2 ). (32) 

47T JO Z 2 JO t 2 



The first term would be easily understood if one notices that the sum of contri- 
butions from real photon emissions and loop diagrams is finite(no infrared and no 
collinear singularity). In this expression the subscript 1(2) designates the electron 
with momentum Pi(P 2 ) and each contribution corresponds to the radiation from 
them. It is clear then by superposing two showers defined with the constraint z > t, 
we get all the contributions in the full phase space. Also the radiated photons are 
parallel to the parent electrons. It is this constraint that allows us to draw a picture 
of jets, a cluster of electrons and photons flowing in the same direction. One should 
notice that in Eq. (|32|) , x B j can be replaced by (1 — Zi)(l — t 2 ). 

Let us apply the above argument to showers with indefinite number of photons. 
Here we will repeat the main points to be argued. In the double cascade scheme, 
parton showers develop with constraint (1 — x) > t, where x is the fraction for the 
electron, not for photons. Then what is the physical quantity that can be calculated 
by the renormalization group equation and what is the combination of the variables 
that equal to this physical quantity? 

For the deep inelastic scattering it is obvious that the answer for the first question 
is the Bjorken variable. Then we study the second question. For this we specify the 
process as 



Here X denotes anything. A virtual electron with pi collides with the virtual photon 
with q and turns into another virtual electron with p 2 . The momentum conservation 
gives 



e(Pi) + 7 (g) -> e(pi) + X + 7 (g) -> e'(p 2 ) + X - X' . 



(33) 



p\ + 2(p 1 q) + q 2 = p\ 



(34) 



In the frame 



P 



0, q+ = -Pi+, qr = 



Pi+ 



(35) 



x b = 



we have 



x Bj = 



2(Piq) 



= x b - 



2(Piq) 




(36) 
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by noting 

2(P ig ) - Xb 2(P iq ) ■ {d7) 

Since |px| = 0(1 — x b )pf, it can be neglected as well as the product of x b and 
p 2 2 /{2{P 1 q)). We find 

x Bj ~ x b (l - yjT^J = x b (l - h). (38) 

In the double cascade scheme p\ > 0. This means that the electron p 2 has a potential 
to emit photons and we must include contributions from showers of the scattered 
electron. 

Next we have to prove that the distribution over x b (l — t 2 ) in fact agree with the 
solution of the renormalization group equation. In the course of the proof one will 
see that this combination of variables is suitable for our discussion. Also after the 
proof, one will find the finite contribution we are looking for. 

The ^-distribution is nothing but the distribution in the single cascade but with 
the constraint (1 — x) >t inside of the integral in exponential. 

D x {n,Q 2 ) = exp J q2 I P(x)(x^ - 1)0(1 -x-t)) . (39) 

The expansion of the distribution over t 2 is 

D t (t 2 ,Q 2 ) =U c (Q 2 ,Ql)5(t 2 ) 

t ? t dxP(x)9(l -x- t)9((l - x)t - Ql/Q 2 )T1 C (Q 2 , tQ 2 )5(t - t 2 ). 

2,71 Jo t Jo 

(40) 

Here U C (Q 2 ,K 2 ) is the probability of non emission with the constraint, i.e. the 
probability in the double cascade scheme. 

II C (Q 2 , K 2 ) = exp (-£- f 1 ^ f 1 dxP(x)9(l -x- t)9((l - x)t - Q 2 /Q 2 )) . 

\ 2ix Jk 2 /q 2 t Jo ) 

(41) 

The moment with respect to t 2 is defined by 

D t (n, Q 2 ) = f 1 dt 2 (l - t 2 )"- 1 A(i 2 , Q 2 )- (42) 
Jo 
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Noting that (1 — t) n 1 ~ 9(l/n — t) for large n, we find 

D t (n,Q 2 ) = exp(^-j Q2 —j o dxP(x)((l - t)- 1 - 1)6(1 - x - t)\ , 

e Q = Ql/Q 2 , t = K 2 /Q 2 . (43) 

We have obtained closed forms for the moment distributions using some approxima- 
tion. Since we know that the product of xt, and 1 — ^2 is the Bjorken variable, the 
moment of XBj, Dd(n,Q 2 ) is the product of the moment D x (n,Q 2 ) with respect to 
Xb and D t (n,Q 2 ) to t 2 , 

D d {n, Q 2 ) = D x (n, Q 2 )D t (n, Q 2 ). (44) 

Next we compare these moment distributions with those obtained by Monte Carlo 
simulation, which will confirm our analysis. First the distribution of Xb will be 
discussed. The analytic form of the moment distribution can be calculated easily. 

D x (n,Q 2 ) = ex V l — j Q2 —j o dxP(x)(x n - 1 - 1)6(1 - x - t)j 

= expf^/ 1 (telog(-)P(l-z)[(l-«)»- 1 -l]y (45) 
\2tiJ^q e Q j 

Figure 2 shows the Monte Carlo results for ^-distribution in the double cascade 

scheme and those obtained from the analytic form Eq. ((45"|) . The agreement justifies 

our formulation. 

Now let us calculate the finite contributions in the double cascade scheme. First 
we notice 

1 . 1 . . 2-t -2 + z 

_ P( 1 _ z) = _ P( i _ t) + __ + __. (46) 

The moments in the scheme are then 
D d (n,Q 2 ) = D x (n,Q 2 )D t (n,Q 2 ) 

= exp(> C- r Q/t dxP(x)(x^-l)) 



2tt Jt Q t Jo 

/a r 1 r 1 fl — t 2 — z\ 
x expf — ^ dt ^ dz [— —) 6(zt - e Q )6(z - t) 

xKl-tr- 1 -!]). (47) 
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The first exponential on the right-hand side is the same as the moments D s (n, Q 2 ) 
appearing in the single cascade scheme. In the latter the integration over z is done 
and 0(€q) terms are neglected. After changing variable t by z, we get 

D d (n,Q 2 ) = D s (n,Q 2 ) 

x exp (|- f Q dz ((-2 + z) \og{z) - 1 + 2 - [(1 - z)^ - 1]) . (48) 
Hence we finally find 

D d (n,Q 2 ) = D(n,Q 2 )D df (n) } (49) 
where the finite term is given by 

D df (n) = ^Hexp^/^z^ + ^logW-A + S-D^l-zr 1 -!]), 

/ a r 1 3 z 

= exp — / dz[\og(z)P(l-z) + ((-2 + z)\og(z)- — + 2--)} 

V Z7T JO 2Z 2 

X[(l-z)»-l-l]). (50) 

This analytic form can be compared with the Monte Carlo results in Fig.2. The 
agreement confirms definitely our discussions. Hence to get the solution of the 
renormalization group equation, D(n,Q 2 ), we have to subtract the finite contribu- 
tion in Eq.(|50|) or calculate the ratio of distributions with different Q 2 . Thus we 
have completed the analysis of the deep-inelastic scattering. 

In e + e~ annihilation, a more complicated situation arises due to its kinematics. 
It is the fact that the squared mass of the virtual photon, q 2 , is calculated by 
momentum fractions Xi,X2 of the electron and the positron. 

q 2 = XiX 2 s, (51) 

where s is the square of the energy in the center-of-mass system. 

In the double cascade scheme we consider the following process allowing both 
electron and positron to radiate photons, 

e-(P x ) + e + (P 2 ) - e-( Pl ) + e + (p 2 ) + X - 7 (g) + X. (52) 
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By this equation we mean that the spacelike virtual electron(pi) and positron^) 
annihilates into the virtual photon(g). In the frame that = (Pi + ,0_,Ot) P 2ll — 
(0 + , P2-1 Ot) a simple algebra gives 

q 2 = ( pi+p2 ) 2 

= pl+pl + 2x bl x b2 (P 1+ P 2 _) + iP \ + Pl Jf 2 ± P f } - 2 Pl • p a , (53) 

2x bl X b2 (P 1+ P 2 -) 

where 

2*1 = -5— , x fe2 = -5—. (54) 

If one neglects the terms of order (1 — x b )t, one finds 

q 2 = x bl x b2 (l-t 1 )(l-t 2 )s, (55) 

where t± —\ p\ \ /s, t 2 —\ pi \ /s. 

In the deep-inelastic scattering one calculates the Bjorken variables in terms of 
the light-cone fraction of the electron in one shower and the squared virtual mass 
of the electron in another independent shower. XBj = x b \{\ — t 2 ). While in the 
annihilation process both of the light-cone fraction and the squared virtual mass 
determined by the generated shower are needed. Hence the moments with respect 
to the variable x b (l —t) are required for one shower. These are given by 

JO t 

X £L f 1 dxP{x)x n ' 1 9{l -x- t)9((l - x)t - Ql/Q 2 ) 

2,71 Jo 

( a ftQ 2 dK' 2 r \ 
x exp J _ J dyP(y)W - 1)6(1 - x - K"/Q*)6((l - x)K" - Qg)J . 

(56) 

By using (1 — t) n ~ 6{l/n — t), it turns out that the moments of Q 2 /s is given by 
the square of D a (n, Q 2 ). 

D a (n,Q 2 ) = f dH n - l D a (i,Q 2 ) 
Jo 

I <y rQ 2 dK 2 r^a \ 
x exp — j( dxP(x)((x(l - f))"" 1 - 1)5(1 -x-t)j, 

e = q 2 /^ (57) 
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The finite contributions D a f(n) is then given by 
D a (n,Q 2 ) = D(n,Q 2 )D af (n,Q 2 ) 

D af (n,Q 2 ) = exp(^- f 1 dx[P{x) hg(l- y/^ -2— log -^^= 
\2tt Jo x 1 + v 1 _ x 

+ 1 T ^log(l-^)](^- 1 -l)), (58) 

for e Q < 1. 



Section 4. The model 

Discussions in the previous sections suggest that it is easy to correct the finite 
contributions in the moment space. For the annihilation process 

D(n,Q 2 )=D a (n,Q 2 )/D af (n). (59) 

Since D a j(n) _1 < 1, its inverse transform D a f(x) is well defined, 

1 r-ioo+c 1 

D af (x) = — / x~ n —--. (60) 
2m J~ioo+c D a f(n) 

We will show how to take into account these corrections in e + e~ annihilation process. 

As a simple analytic form is not easy to get for D a f(x), we make the following 

approximation 

D af (x) ~ D A af {x) = C5(l - x )+e(x f - x)F(x), (61) 
where the function F(x) is given by 



P(x) ln(l - y/x) - 2^-^ In v "'' ""' 



x 1 + Vl — x 

+ ^^ln(l-a;)lexpf-^ln 2 (l- a ;)V (62) 
2 J V 2-7T / 



Here 1 — Xf <^ 1 is assumed and the constant C is fixed by the requirement that 

dxD^ f {x) = 1. (63) 



o 



The parameter Xf is fixed by the condition that a photon with energy fraction 
I — x f cannot be observed in the detectors. In usual experiments the limit of the 
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measurable energy is in the order of 10 MeV. From this fact in the program we set 
xj = 1 — 1CT 4 . If 1 — Xf is extremely small, the approximation of D a Jx) for D a f(x) 
is not justified. Then we add the following algorithm before developing showers. 

At Ql, we choose x according to D^(x). If x ^ 1, a photon with the momentum 
fraction 1 — x is emitted parallel to the electron and the energy fraction of the 
electron is x. In the generator Q 2 S S> Ql is assumed. 

Next we will present results by the improved QEDPS and compare them with 
the analytic results. Since the coupling constant a is small, we make a comparison 
on the integrated radiator over small 1 — x, that is, 



In the Monte Carlo we count the fraction of events with s e ff ec u ve > xs. In the 
simulation we assumed the following values for the parameters. 



The results are shown in figure 3. 

In analyzing the experimental data people use very sophisticated radiators, which 
involve the higher order corrections and others. A review on these will be found in 
ref.0. We call the radiator given there Rs upp i{x, Q 2 )- This and R(x, Q 2 ) are plotted 
in figure 3, in which one still finds a small discrepancy on the order of 1%. We should 
note, however, that 




(64) 



(65) 



a = 1/137, Q 2 S = l(T 6 GeV 2 ,Q 2 = 10 4 GeV 2 ,Q 2 = HT 8 GeV 2 . 



(66) 




(67) 



while 




(68) 



We introduce a modified Rm(x, Q 2 ) as 



R M (x,Q 2 )=R(x,Q 2 )C a . 



(69) 
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Then we plot Rm(x, Q 2 ) and R,s U ppi{x, Q 2 )C a to find a satisfactory agreement whose 
accuracy is less than 0.1%. Hence the normalization is fixed by C a instead of unity 
in the shower. 

Section 5. Summary and discussions 

In this paper we present the improved version of QEDPS, a generator for radia- 
tive corrections to the processes in e + e~ annihilation. This new model has been 
obtained from the old one, that has been proposed by the authors a few years ago, 
by correcting the finite contributions originated from our shower algorithm. It was 
shown that these contributions appear when the infra-red cutoff is Independent. 

In the present work we did not take account of the running effect for the QED 
coupling a. If one applies this study to improve the finite contributions in QCD 
showers, the running effects should be considered. We would like to discuss this 
point in a separated paper. 

Acknowledgements 

We would like to thank our colleagues of KEK working group (Minami-Tateya) 
and in LAPP for their interest and discussions. This work has been done under 
the collaboration between KEK and LAPP supported by Monbusho, Japan(No. 
07044097) and CNRS/IN2P3, France. 



17 



References 

[1] G. Bonvicini and L. Trentadue, Nucl. Phys. B323(1989)253, 
S.Jadach and B.F.Ward, Comput. Phys. Commun. 56(1990)351, 
Darmstadt-Siegen Collaboration (H.D. Dahmen et al.), Z.Phys. C50 (1991)75, 
S. Jadach, E. Richter-Was, B.F.L. Ward and Z. Was, Comput. Phys. Commun. 
70(1992)305, 

M. Caffo, H. Czyz and E.Remiddi, II Nuovo Cim. A105(1992)277; Int.J.Mod. 
Phys. C4,No. 3(1993)591. 

[2] J. Fujimoto, Y. Shimizu and T. Munehisa, Prog. Theor. Phys. 90(1993)177. 

[3] J. Fujimoto, Y. Shimizu and T. Munehisa, Prog. Theor. Phys. 91(1994)333. 

[4] R. Odorico, Nucl. Phys. B172(1980)157, 

G. Marchesini and B.R. Webber, Nucl. Phys. B238(1984)l. 

[5] K. Kato and T. Munehisa, Phys. Rev. D36(1987)61, 

K. Kato and T. Munehisa, Comput. Phys. Commun. 64(1991)67, 
T. Munehisa and H. Tanaka, Int. J. Mod. Phys. A6(1991)989, 
K. Kato, T. Munehisa and H. Tanaka, Z. Phys. C54(1992)397. 

[6] K. Kato and T. Munehisa, Phys. Rev. D39(1989)156. 

[7] V.N. Gribov and L.N. Lipatov, Sov.J.Nucl.Phys.l5(1972)438, 

E.A. Kuraev and V.S. Fadin, Sov.J.Nucl.Phys.41( 1985), 466; DESY Internal 
report L-Trans-297,1985, 

J.P Alexander, G. Bonvicini, PS. Drell and R. Frey, Phys. Rev. D37(1988),56, 
O. Nicrosini and L. Trentadue, Phys. Lett. 231B(1989)487, 
J.Fujimoto et al., Prog.Theor.Phys.Suppl. No.l00(1990)l. 



18 



Figure Captions 

Fig.l Moments of the structure functions. The solid curve is the correct re- 
sults obtained by the renormalization group equation. The dashed curve is given 
by the analytic form of D s (n,Q 2 ). Diamond(circle) data points are Monte Carlo 
results in the single cascade scheme for Q 2 = Ql(Q 2 s 3> Ql). Parameters are 
Q 2 = 10 4 GeV 2 ,Q 2 = io-&GeV 2 ,a = 1/20. For the Monte Carlo data of Q 2 > Q 2 , 
we choose Q 2 = l(T 8 GeV 2 . 

Fig.2 Moments of the structure functions in the double cascade scheme. The 
solid(dashed) curve is calculated analytically by D x {n,Q 2 ) (D x (n,Q 2 )D t (n,Q 2 )). 
Diamond(circle) data points of the Monte Carlo simulation corresponds to x b (x B j = 
Xb(l — t)). The long-dashed curve presents the results of D(n,Q 2 ). Parameters are 
the same as Fig.l. 

Fig. 3 Differences between the Monte Carlo data and analytic expression for the 
integrated radiator, J* dyR(y, Q 2 ). Parameters are the same as Fig.2 except a. Here 
a = 1/137. One million events are generated. 
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